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Abstract 


We discuss specific, recent advances in the analysis of an experiment to test the 
Equivalence Principle (EP) in free fall. A differential accelerometer detector with two 
proof masses of different materials free falls inside an evacuated capsule previously 
released from a stratospheric balloon. The detector spins slowly about its horizontal axis 
during the fall. An EP violation signal (if present) will manifest itself at the rotational 
frequency of the detector. The detector operates in a quiet environment as it slowly 
moves with respect to the co-moving capsule. There are, however, gravitational and 
dynamical noise contributions that need to be evaluated in order to define key 
requirements for this experiment. Specifically, higher-order mass moments of the 
capsule contribute errors to the differential acceleration output with components at the 
spin frequency which need to be minimized. The dynamics of the free falling detector (in 
its present design) has been simulated in order to estimate the tolerable errors at release 
which, in turn, define the release mechanism requirements. Moreover, the study of the 
higher-order mass moments for a worst-case position of the detector package relative to 
the cryostat has led to the definition of requirements on the shape and size of the proof 
masses. 


Introduction 


The scientific goal of the experiment is to test the equality of gravitational and inertial 
mass (i.e., to test the Principle of Equivalence) by measuring the independence of the rate 
of fall of bodies from their mass and composition. The experiment is accomplished by 
measuring the relative displacement (or equivalently acceleration) of two falling bodies 
of different materials which are the proof masses of a differential accelerometer. The 
experiment goal is to measure the Eotvos ratio 6g/g (differential acceleration/common 
acceleration) with an accuracy about two orders of magnitude better than presently 
achieved 111 . The experiment is a null experiment in which a relative displacement 
different from zero, at the spin frequency, between the proof masses will indicate a 
violation of the Equivalence Principle. 
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In summary, the experiment consists in taking differential acceleration measurements 
with a high- sensitivity detector (the sensor) during free fall conditions lasting up to 28 s 
in a disturbance-free acceleration environment. A capsule is first released from the 
balloon at an altitude of typically 40 km and the detector is released from the top of the 
capsule immediately afterwards. During the measurement phase, the sensor free falls 
inside a 2-meters-long (in the vertical direction) evacuated cryostat (contained inside the 
capsule) that is falling simultaneously in the rarefied atmosphere 111 . 

By falling in vacuum inside a co-moving capsule, the noise acceleration level can be 
kept to a negligible level while the signal strength in free fall, i.e., the full-strength 
Earth’s gravity, is increased by 3 orders of magnitude with respect to the signals available 
to EP experiments conducted on the ground. The free fall technique, therefore, combines 
some of the advantages of the space-based tests with the accessibility and reusability of 
ground experiments. 



Detector before 
release 


Spin system 


Sensitive axis 
(rotates with detector) 

Cryostat 


Figure 1 Schematic of capsule with detector attached (before release) to the spin up 
system 

The detector has two sensing masses of different materials with their centers of mass 
(CM) as coincident as possible in order to reduce the effect of gravity gradients. The two 
masses are constrained in a way that they can move at low frequency along one axis (i.e., 
the sensitive axis that is perpendicular to the symmetry axis of the sensor) while they are 
rigidly constrained about the other two axes. The displacement of each sensing mass 
along the sensitive axis is detected by capacitive pickups which are parts of a 
measurement bridge. 
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One key feature is that the detector is spun before release about the symmetry axis 
that lies on the horizontal plane. The rotation provides the modulation (at the spin 
frequency) of any possible violation signal. Moreover, the spin provides gyroscopic 
stabilization of the detector during the fall and separates the signal from key noise 
components generated by the detector dynamics. 

Frequency separation between the signal and the largest noise components is essential 
to the success of the experiment. The noise components that, unavoidably, appear at the 
signal frequency need to have a strength well below the detector threshold sensitivity. In 
the following, we analyze the detector dynamics to identify the frequencies of noise 
components that can be separated from the signal thanks to the appropriate selection of 
the inertia characteristics of the detector package. We also study the higher-order mass 
moments generated by the capsule mass and acting on the sensing masses. The higher- 
order mass moments contain some components at the signal frequency and care must be 
taken in designing the sensing masses in order to reduce the strength of those components 
to a negligible level. 


Analysis of Higher-order Mass Moments 


Overview 

The sensing mass (test body), falling inside the capsule, is subjected to non-negligible 
gravitational attraction by the capsule. The higher-order gravitational potential plays a 
key role because of the elimination of the zero-order potential due to the Equivalence 
Principle. The model of perturbing gravitational mass consists of a spinning test body 
inside a capsule that is a hollow cylinder covered with flat caps. The test body is released 
at the symmetry axis of the cylinder, and deviates from the axis during its fall. Our goal 
is to compute the gravitational force and torque acting in the neighborhood of the fall 
trajectory. The fact that both the test mass and the capsule are closed finite bodies, 
increases the complexity of the problem. We are mostly concerned here with the 
harmonics of the force/torque at the modulation frequency as seen in the reference frame 
of the detector. In other words, our model should evaluate the force/torque in the rotating 
body frame of the sensing mass. For that purpose we built a semi-analytical model that 
can handle any configuration of test mass as well as any additional mass attached to the 
distributed mass of the capsule. Moreover, we exploited the fact that the size of the test 
mass is smaller than the cylinder radius to derive an asymptotic analytical solution. 

There are at least three ways to approach the computation of the gravitational 
attraction between the capsule and the test body. The first approach is to compute the 
force between each capsule mass element and a test body mass element, and to perform a 
double summation on these forces. If N is the number of capsule mass elements, and 
N B is the number of test body mass elements, then the cardinality of the computation 

(that drives the computational load) is N ®N B . As it will become apparent later on, the 

mass discretization resolution required to the test body is very demanding and, 
consequently, the cardinality becomes prohibitively high with this approach. The second 
approach is a double integration over the bodies. The drawback of the first approach is 
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the heavy computational effort, especially if the computation is needed as an online 
computation in a dynamical process. The weakness of the second approach is in the 
cumbersome analytical computation. The integration will require an asymptotic series 
expansion of complicated functions and asymptotic series introduce a truncation error 
into the computation. We adopted a third approach for the gravitational model. 

We consider the capsule as a discrete ensemble of lumped mass, where the resultant 
force and torque acting on the test body are the cumulative force and torque due to the 
summation of all the capsule mass elements. The gravitational potential of the test body 
is expanded according to its increasing-order inertia characteristics. In other words, the 
interaction is between a finite body (the test mass) and a point mass. This approach is the 
result of a tradeoff between the previous approaches and its cardinality is N . The main 
advantage is the flexibility of modeling any capsule shape and the strongly reduced 
computational load. 

The purpose of the following computation is to derive a closed-form analytical 
expressions, for the force and torque acting on the spinning test mass due to a capsule 
mass element. 

Gravitational Model 

The gravitational potential between finite size bodies, is: 

dM 

V = -G I f ~dM (1) 

m m r ° 

o B 

where 9^ is the test body (proof mass), and d( 0 represents the attracting bodies, namely, 
the capsule or any other perturbing body. 

For the purpose of simplicity we will proceed with a formulation of the test body and 
a single element of the attracting mass, M i . 

The gravitational potential at a representative element mass M i due to the test body 

dM g _ 

is: V(M .) = -GM . f , where r is the radius vector between a mass element 

X r 

of the test body and dM a as shown in Fig. 2. 
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Figure 2 Gravitation model for test mass and outer attracting bodies 

Assuming a sufficient discretization, the total force acting on the test body due to all 

N 

external mass elements is F B & where N is the number of external mass 

1 = 1 

elements. We first erect the following reference frames. The capsule frame, denoted by 
{X,Y,Zj is attached to the capsule. X,Y lie on the equatorial plane and have an arbitrary 
azimuth while Z coincides with the symmetry axis of the capsule/cylinder. The test body 
frame, denoted by {x,y,zj , is attached to the test mass. The gravitational potential will be 
expressed in body frame. The inertia coefficients of the body are constant in this frame. 

The derivation consists of the following steps: (1) substituting the inverse radius 
approximation into the potential; and (2) carrying out the integration. The position 
vectors in body frame are as follows: 

R = (le x + me v + ne z )R and p - xe x + ye v + ze z 

where {e x , e v , e z } are body unit vectors, and {/, m, n) are direction cosines between R and 
the body axes. 

The resulting potential is a summation over all mass elements M; (only the first 4 
terms are shown for the sake of brevity): 
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N GM x M 2 GM x M 2 

V = — — V - 4 \lx + my + nz\ 

i = 1 R R z 


GM 


2R 


GMi i i i 

^-[Z(5/ —3)Jxxx + m (5 m ~3) Jyyy +n(5n —3)J 

27? 

+ 3m(5l 2 -l)J xxy + 3n(5l 2 -1) J xxz + 3l(5m 2 -1) J xyy 
+ 3 n (5m 2 - 1) J yyz + 31 (5 n 2 - 1) J + 3m (5 n 2 - 1) J yzz 
+ 6lmnJ xvz ] + 0(^)} 


±[(3l 2 -1)J XX + (3m 2 -l)J yy + (3n 2 -l)J zz + 


6(lm J xy + InJ xz + mn J yz)\ 


'xyz 


7?~ 


( 2 ) 


where j xPy q z r are the inertia integrals defined as: 


'x p y q z r 


- f x p y q z r dm 


Mi 


(3) 


Expressions for the forces (shown in the following) are obtained by derivation of the 
potential. Torques (not shown here) are obtained by integrating over the proof mass the 
torque acting on a mass element of the body. 


F x = Gl M,.{^ + ^%[(3/. 2 - l)x + 3/-m • y + 31- n- z 
i = 1 Rf Rf 

+ ' 4 K 5 ^ ~hy ~ l zz + + 3m.(5/3 -l)I xy + 3«.(5/ ; 2 -1)1^ + 15/^1*] 

R i 

+ ^5 [ 2 (35/ / 4 “ 30/ / 2 + 3 ) J x« + 2 (7m / 2 " 3)l i m i J m + \ {ln i ~ 3) h n i 
+ y (7 Z,. 2 - 3 )l i m i J xxy + y (7/3 - 3) /• n. J xxz + 1 (35/ 2 m3 -5/3 - 5m- 2 + 1) J xyy 
+ y (7m3 -1)/^. J yyz + |(35/. 2 „. 2 -5/3 -5n. 2 + 1)J XZZ + ^(7/». 2 -l)/. m .J yzz 


+ 3(7 1- 1 - 1 ) m-n- J xyz J + H D.T\ 
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N M B m . T ? 1 

= G 2 * ? I 31m- x + (3m- - l)y + 3m -n- z 

i = 1 Rr*'- 1 1 J 

15 2 

+ ^ [(— - I xx - + 21 yy )m^ + 3 (5 - 1) I X y +15 Ixz 

A ’/ 


^xz + 3n i (5m i 2 -l)I yz ] 




+ tf 5 2 


f (35m3 - 30m 2 + 3) J yyy + |(7n3 - 3 )m ; y J zzz 


1^5 1 2 m 2 -5m 2 -5l 2 + l)J xxy + ^ (7/4 - l)m ; n,.J xr - + y(7m- 2 - 3)/.#».J xyy 

- y(7m ; 2 _ y > m i n i J yyz + y( 7n ; 2 -^h^xzz + |(35m3 W/ 2 -5m- 2 - 5/7 - 2 + l)J yzz 

- 3(7 m 3 - 1) / f n- J xyz ] + H .O.T} 


= G 2 M i { ^ 3/ n x + 3m -n- y + 3(n 2 - 1) z 

i-i y R i 

+ ' 4 [(J T. - I xv - I vv + 2I s )n. + 1 5l-m-n- I xy + 3/-(5n 2 -1)I W 

y- z 


1 5 2 5 2 1 

+ — -3)/- n | J XC( + -(7m- - 3)m- n- + -(35n- 1 - 3Un-“ + 3)J SZ 
+ y (7 1? - ] ) m i"i } xxy + 1(35 l?n 2 - 51 2 - 5 n 2 + 1)J W + y (7m 2 - 1) /• n- 


n 2 -l)\ xz + 3mpn 2 -l)\ yz \ 
^(35n- 4 -30 n 2 + 3);y z 


|( 7 /. 2 -Dm.n-J^ + f (35/-V -5/ 2 -5, 
■ |(35/n. 2 n . 2 - 5m. 2 - 5„ ; 2 + 1) + y (7y 


y ( 7 «/ 2 - HZ, «/ A z + y Vn 2 - 3)m.n. J yzz 


+ 3(7 n} - 1) l-m- J rv ,] + H.O.T} 
where H.O.T stands for higher-order terms and 

2 2 : 
r. = (I yy + I ZZ - 2 Ax A + 4 XX + Az - 2 AA m j + (I XX + At - 2 AzA/ 

Note that the force has the following order of magnitude: 

2 3 

F oc 2G Mi ^ B {0(l) + 0(-A) + 0(^y-) + 0(^\) + H.O.T} 

i Ri R t R i 2 R 2 

Where L B is a characteristic length of the test body with — « 1 , and — « 1 . 

R , R i 


2 ' 
2 l yy) m i + (Ax + At — 2 Az)y 
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Minimizing the higher-order effects 

The goal is to minimize the gravitational forces acting on the test mass. The force 
equations reveal that the dominant term is the term corresponding to the second-order 
inertia. The offset term (dipole) is equal to zero for a reference frame centered at the 
body CM (see also later on). The direct way to minimize the forces is to require equal 
second-order moments of inertia. Then, the residual forces are due to mass imperfections 
and to higher-order (and smaller strength) inertia integrals. The purpose of the following 
discussion is to explore the effect of each term for a spinning test body. 


First, we evaluate the mixed-inertia terms due to an imperfection. We assume that the 
imperfection is due to a disturbance in mass and/or length. The nominal test body could 
be a perfect cylinder, with radius R B and length L B . In order to minimize the force we 
require that I xx = I VY = / zz . This constraint translates into a given aspect ratio of the 

cylinder as follows: L B =-^3 R B . Given the mass density, the mass or the size of the 


cylinder are now functions of a single free parameter. For example, given the mass, M B , 

l 


/u/ \ 


and the density, p B , the length is: L B = 


3 M 


B 


\Jlp B j 


The similarity dimension of the k-order inertia is [/] = 



. Thus, the perturbed k- 


order inertia is 51 = i} 5M + kMl} 1 5L , or 


'81 


8M 

+ k 

8L 

— 

= 





_ M _ 


L 


Next, we consider the effect of the spin. The test body is spinning about its x-axis 
with a frequency that will be regarded as a fundamental frequency, or PI (period one). 
We are primarily concerned with PI because it is the frequency of the measured signal. 
For this purpose, we will analyze F y by substituting the direction cosines, shown below, 
into the force expression. 

The periodicity in time is introduced through the direction cosines that represent the 
orientation of the body frame relative to the capsule frame. The transformation between 
the two is: 


'x 

y 

— 

'I 0 
o cd 

0 ■ 

sd 

' X ' 
Y 

(5) 

z 


0 

1 

< 3 > 

cd 

z 



The resulting direction cosines are: 
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x X 

l = — = — ( independent of B ) 

R R 

y YcB+ZsB 

m = — = ( 6 ) 

R R 

_ z _-Y sB+ ZcB 
H ~ R~ R 


l becomes a cyclic function only if the body deviates from the X-axis. 



Figure 3 Reference frames 

We compute the forces and torques acting on an imperfect proof mass placed at the 
worst location expected inside the capsule for conservative value of the wind shear (that 
is a point at -0.5 m below the capsule equator and 0.1 m off the centerline. The capsule 
consists of a hollow cylinder, covered by two flat caps with a total mass of 500 kg. The 
test body mass is 1 kg with sizes: R B = 0.0412m and L B = 0.0713m. The geometrical 
and mass inhomogeneity errors, based on realistic assumptions, are: 6L/L = 10' 4 and 
6M/M = 10‘ 4 . 

In the following plots (Figs. 4 and 5) we also show non-zero dipole terms (first-order 
inertia) to represent the fact the CM of one proof mass does not coincide with the CM of 
the other, thereby originating a non-null (gravity gradient) term when the difference of 
accelerations is taken between the two proof masses. 
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Figure 4 Ordered forces and spectra at expected worst location within the capsule 
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Figure 5 Ordered torques and spectra at expected worst location within the capsule 
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Similarity analysis 

Another important issue is the role played by the test body mass and size. For this 
purpose we perform a similarity analysis, as shown in the following. We have formulated 

the gravitational potential as an asymptotic series in 1/R n . Each term in the series has an 
inertial tensor of degree n - 1. For example, the first term is of order n = 1 and its inertia 
tensor is of order 0 (the test body mass). The third term is of degree n = 3 and its inertia 
tensor is of order 2. This non-uniformity, however, introduces a problem when one 
computes the forces for a particular test mass and wants to deduct the forces for a 
different size of test mass. Mathematically speaking, if F(M m ) and F(M B2 ) are the forces 
due to two different test masses, then the ratio between the forces is not a homogeneous 
function, that is, F(M BI )/F(M B2 ) * g((M BI /M B2 ) m ) where g is a function and m is the degree 
of the homogeneity. Note however that each term of the series is homogeneous by itself. 
For example, the first term is homogeneous of degree one, that is, the ratio between 
forces equals the ratio between the masses. 

We will distinguish between two situations. The first situation is when the difference 
in masses is due to a different density. Since the forces are homogeneous of degree one in 
density, the ratio between the forces is equal to the ratio between the masses and the 
acceleration is constant. In a more interesting situation the density is the same, and the 
different mass is due to different sizes. In this case we need to apply a different similarity 
to each term. Let L and M be scaling factors for the length and the mass, respectively. If 

the density of different test bodies is the same, then LocM . Let m = n - 1 be the order 
of the moment of inertia, then its similarity dimension is L m M , that is, M 1+m/3 or L m+3 . 

The acceleration similarity is M m/3 or L m . The following table summarizes the 
similarity dimension for each term in the potential. 


Table 1 Similarity relations for the gravitational potential term 


Potential Order 

Inertia Order 

Inertia 

Acceleration 

n; 1/R n 

m = n - 1 

Similarity 

Similarity 

1 

0 

M (L 3 ) 

1 ( 1 ) 

2 

1 

M 4 ' 3 (L 4 ) 

M" 3 (L) 

3 

2 

M 513 (L 5 ) 

M 2,3 (L 2 ) 

4 

3 

M 2 (L 6 ) 

M (L 3 ) 

5 

4 

M m (L 7 ) 

M 4 ’ 3 (L 4 ) 


To summarize, given the force (or the acceleration) on a particular test body, we can 
deduct the force (acceleration) on another geometrically-scaled test mass. The 
acceleration ratio (of two different test masses with the same density) versus similarity 
dimensions (i.e., mass and length) is illustrated in Fig. 5. 
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Figure 6 Scaling of test mass acceleration for masses with same density 

The maximum potential order considered in our analysis is the hexadecapole (n = 5) 
which involves the fourth-order inertia integrals. We computed numerically the forces 
associated with the hexadecapole for a body with equal second-order inertia integrals (in 
order to minimize the quadrupole term). The results show that, at the worst expected 
location inside the capsule, the maximum acceleration associated with the hexadecapole 
term for a 1-kg proof mass is less than 10' 16 g. In other words, for the accuracy goal of 
this experiment there is no need for belted cylinders iv (which would reduce the 
hexadecapole component even further). Proof masses with equal second-order inertia 
integrals (or alternatively moments of inertia), construction accuracy of order ten 
microns, and sizes smaller than about 10 cm are sufficient to make the contributions of all 
the higher-order gravity terms negligible. 

Detector Dynamics 

The instrument package dynamics also affects the differential accelerometer output. 
We derive the differential equations of motion for the detector in free fall (without 
gravity perturbations for the time being) and carry out numerical integrations for cases of 
interest. Those cases are associated with non-null offset positions of the centers of mass 
of the proof masses with respect to the CM of the instrument package and initial 
rotational velocity errors orthogonal to the spin velocity in order to evaluate the influence 
of those parameters on the differential output of the accelerometer. 
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Equations of motion 

The sensor in its current, preliminary design (see Fig. 7) is composed of two sensing 
masses (A and B) having the CMs as coincident as possible with the CM of the external 
case. One proof mass has a dumbbell shape, while the other is a hollow cylinder. These 
two masses are constrained to the case C by means of elastic springs and they can pivot 
about an axis that is parallel to the axis of symmetry of the sensor. The whole detector is 
spun about the symmetry axis x. 


C 




A 


B 

X 

r 

1 — \ 

|L- 

- 

I 1 1 

w 

T” 


spin axis 


pivot axis 


Figure 7 Schematic of three-body sensor 

The motion of the sensor is essentially a rotation of the sensing masses about the 
pivot axis. The translational motion of the sensing masses, which depends on the high 
lateral stiffness of the springs is much smaller than the rotational motion. A rotation of 
the proof mass causes a variation of the distance between the sensing mass and the 
external case that changes the output capacitance. 

Each body is modeled with six degrees of freedom, i.e., three translations and three 
rotations defined by Euler’s angles. Elastic forces and torques are applied to each proof 
mass at the attachment points of the constraining springs. The transformation matrices 
have been written for each body in order to develop the equations of motion for the 
different bodies. The body coordinate systems are centered in the body’s geometrical 
center, and are fixed with each body. The position of the CM is defined by three 
coordinates (x CMA , y CMA , z CMA for body A) in the body reference frame. Each body frame 
has the x axis along the axis of symmetry which is also the spin axis. The y and z axes are 
radial axes that rotate with the body. The y axis is the sensitive axis of the accelerometer. 

Translational accelerations 

In the following we describe the procedure for deriving the equations of motion of 
one proof mass (A). The extension to a second proof mass is straight forward v . The CM 
vector in each body’s reference frame is: 
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Body A: {r A } = (x CMA y CMA 

z cma} 

Body C: {r c } = (x CMC y CMC 

Z CMcl 

The translational acceleration expressed in the inertial reference is obtained by using 
the formula: 

K}-[Mh> +2 [Mfe> + Ko]-fc}+{«S} (7) 

where {a A } is the acceleration of body A in the inertial frame, [R A0 ] is the rotation matrix 
from A to inertial frame, {r A } is the coordinates vector of body A CM in the A frame, 
{R a 0 } is the translation vector expressed in inertial coordinates, and ( ) indicates matrix 
multiplication. Same-structure equations are used for the additional bodies. 

Elastic Forces 

To evaluate the elastic forces we define first the points where the springs are attached to 
each body. The position of these points are expressed in body frame. The distance 
between connected points is then computed from the actual position during the motion, 
projected in the A frame, and multiplied by the stiffness vector {k xA , k yA , k zA }. By 
following this procedure, it is possible to assign different stiffness to each degree of 
freedom. The elastic forces are then projected in the inertial coordinate system as follows 


{F 1A } = [R ao ] • ([K] . ({p A1 } - [T oa ].{T co . Pci } - {1,})) 


{F 2A } = [Rao] • ([K] . ({p A2 } - [T 0A ].{T C0 .p C2 } - {1 2 })) (8) 

{F A }={F 1A } + {F 2A } 

where [K] is the spring stiffness matrix, {p A } and {p c } are the vectors that define the 
positions of the spring attachment points, and {l,} and {1 2 } are the natural lengths of the 
springs. The operator “T 0A .x” indicates the combination of a multiplication by the 
rotation matrix [R 0A ] applied to the vector {x} plus the translation of the vector {R A A } 
which yields: 

{T oa .x} = [R oa ].{x} + {R a a } 

The expression T 0A .(T C0 .p cl ) projects the coordinates of the point p ci (expressed in the 
body C frame) onto the body A coordinate system. Subsequently, the force is projected 
onto the inertial frame using the [R A0 ] matrix. A similar procedure is used to evaluate the 
forces acting on any other body pairs. 

The equations of translational motion for the two bodies A and C finally yield : 


m A {a A } - {F a } = 0 


me {a c } - {F c } = 0 


( 9 ) 


14 

15 



2003 NASA/JPL Workshop on Fundamental Physics in Space, April 14-16, 2003, Oxnard, California 


where m A and m c are the masses of body A and C, a A and a,- the accelerations, and F A and 
F c the elastic forces: 

{F a } = {F 1A } + {F 2A } 

{F c } = {F 1C } + {F 2C } 

The subscripts 1 and 2 identify the two springs that connect body A to C. 

Rotational accelerations: 

The angular velocity of each body is computed by using the rotation matrices that 
transform the coordinate system from the inertial to the body frame and, conversely, 
through the opposite transformation. The rotational velocity matrices of bodies A and C 
are derived by using the Cartan’s formula as follows: 

khm-M 

r • i (10) 

KH*ocH*co] 


where [R 0A ] and [R A0 ] are the rotation matrices from the inertial coordinate system 
(denoted by 0) to the body reference frame of A and vice-versa. [R oc ] and [R co ] are the 
correspondent matrices for body C. It should be noted that the former expression lead to 
the skew symmetric matrix of the angular velocity from which the components of the 
angular velocity vector {oo A } can be readily extracted. 

Elastic torques 

The elastic torques acting on each body are computed by using the expressions for the 
locations of the attachment points and elastic forces previously defined. The expressions 
of the torques in the respective body reference frames are as follows: 


{T a } = {p A1 } x ([R oa ] . {F 1A }) + {p A2 } x ([R oa ] . {F 2A }) 


{Tc} = {Pcil X ([Roc] -{Fi C }) + {PC2} X ([Roc] • {F 2 c}) 

where x indicates the external product of vectors. 

Invoking Euler’s equations, the rotational equations of motion yield: 


( 11 ) 


[1a]-{"a} + ["a]-[1a]'{"a}-{Fa} = 0 

[/c]-K} + K]-[/c]-K}-{f c } = o 


where {oo A } and [oo c ] represent the angular velocity vector and matrix, respectively, for 
body A (and similarly for body C). [I A ] and [I c ] are the inertia matrices for body A and C 
which, assuming principal axes, have diagonal forms. The equations of motion for the 
complete three-body detector are evaluated in the same way as for the two-body example 
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used previously leading to 3 vector equations for the rotational dynamics and 3 equations 
for the translational dynamics. 

Numerical cases 


Figure 8 shows the configuration of the 3-body sensor utilized for the numerical cases 
and the positions of the attachment points for the three-body sensor. 



Figure 8 Attachment points of the proof masses 

The points p CA1 , p CA2 , p CB i and p CB2 are the points of the external body C connected to 
the points p A1 , p A2 , p B1 , p B2 , respectively. In this detector configuration, four springs are 
used to connect the two bodies to the external case (i.e., two springs per proof mass). 
The numerical value adopted for the numerical case shown in the following are as 
follows: 

m A = 1 kg, m B = 1 kg, m c = 30 kg; 

I Ax = I Ay = Iaz = 0.17 kg-m 2 ; I Bx = I By = I Bz = 0.17 kg-m 2 ; 

I Cx = 0.95 kg-m 2 , I Cy = I Cz = 1.46 kg-m 2 ; 
k Alx = 45000 N/m, k Aly = k Alz = 35000 N/m; 
k A 2 x = 45000 N/m, k A2y = k A2z = 35000 N/m; 
k A6x = 61.68 Nm/rad, k A6y = k A6z = 1000 Nm/rad; 
k Blx = 45000 N/m, k Bly = k Blz = 35000 N/m; 
k B2x = 45000 N/m, k B2y = k B2z = 35000 N/m; 
k B6x = 61.68 Nm/rad, k B6y = k B6z = 1000 Nm/rad. 

The initial conditions at release are: spin rate co x = 1.885 rad/s (0.3 Hz), angular error 
at release oo y = 10' 2 rad/s; and initial nutation angle = 0. The CM offset errors are (with 
the reference frame placed at the geometrical center of body C): 

x C ma = 10' 6 m, y CMA = 0, z CMA = 10' 6 m; 

X CMB = 0; YCMB = 0; Z CMB = 0; 

X CMC = O’ y C MC = 0’ Z CMC = 0 
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From the point of view of the rigid-body dynamics, a detector with the above inertia 
characteristics is a minor-inertia-axis spinner that exhibits a prograde precession. 

The natural elastic frequency of the detector are shown in Figure 9. The lowest 
frequency (at 2.999 Hz) corresponds to the differential torsional mode of the proof 
masses while the next one (at 3.486 Hz) corresponds to the common-mode torsional 
frequency. 



2.999 

3.486 

4.349 

4.399 

26.779 

27.217 


42.108' 

42.554 

43.489 

44.075 

48.344 

50.062 


Figure 9 Natural frequencies 
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-2x10 

-4x10 
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Figure 10 y A ' Yb ( m ) i n body frame vs. time (s) 
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Figure 11 FFT of y A - y B (m) vs. frequency (Hz) 



Figure 12 Zoom of y A - y B (m) FFT vs. frequency (Hz) 

Figure 10 is a time profile of the differential displacement between the two proof 
masses, that is measured by the capacitive pickups. The displacement is then converted 
into acceleration through the accelerometer’s transfer function. The initial conditions 
adopted for release are conservatively large and no post-release phase with oscillation 
damping was included in this simulation leading to conservative values of the differential 
displacements. Here we are mostly concerned with the frequency content of the 
differential displacement which is shown in Fig. 1 1 and, with an expanded view, in Fig. 
12. Figure 11 shows the separation between the oscillations related to the rigid-body 
dynamics of the instrument package (at low frequency) and the first elastic frequency of 
the accelerometer. More importantly. Fig. 12 shows that the precession frequency can be 
chosen so that the harmonics related to the rigid-body dynamics do not overlap with the 
spin frequency (which is also the signal frequency). 
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The precession of the instrument package will be caused by an angular rate error 
(perpendicular to the spin axis) at release or an imperfect inertial balancing of the 
instrument package. Effects of the precession are detected by the accelerometer if the 
centers of mass of the proof masses do not coincide with the center of mass of the 
instrument package. If we choose the moments of inertia of the instrument package so 
that the precession frequency is non-commensurate with the spin frequency, then there 
are no dynamics-related harmonics at the signal frequency. Consequently, we must 
simply guarantee that the precession-related accelerations are smaller than the end-of 
scale of the detector. These conditions are met for a prograde precession (minor-axis 
spinner), angular rate errors at release smaller than about 0.1 deg/s, and realistic 
construction accuracies of the detector. 


Conclusions 


Our analysis concludes that the gravitational perturbations, acting on the test masses 
and due to the capsule gravity, can be reduced to within the limit required by the 
experiment if the test mass design abides to simple rules as follows: 

a) The test masses must be smaller than a characteristic size of about 10 cm; 

b) The second order principal moments of inertia must be equal within construction 
tolerances (6L/L < 10' 4 ); 

c) The density uniformity of the test mass must be within 0.01% (6M/M < 10' 4 ). 

There is no need for belted cylinders (used in the STEP satellite experiment) for the 
accuracy goal of our experiment. 

With reference to the detector free-fall dynamics, the ratio of the moments of inertia 
of the instrument package must be such that the body precession frequency is non- 
commensurate with the spin frequency. Similarly, none of the elastic frequencies of the 
detector must overlap with the spin frequency. These conditions guarantee that no 
acceleration components related to the detector elastic and rigid-body dynamics overlap 
with the signal frequency. 
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